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The Decay of Surface Waves produced by a Superposed Layer 

of Viscous Fluid. 

By W. J. Hakrison, B.A., Fellow of Clare College, Cambridge. 

(Communicated by Prof. J. Larmor, Sec. E.S. Eeceived May 13, — Eead 

June 17, 1909.) 

§ 1. In a paper published recently by Prof. E. W. Wood, on " The Mercury 
Paraboloid as a Eeflecting Telescope,"* he notes the complete damping of 
ripples in the mercury surface produced by a thin superposed layer of 
glycerine. It has been suggested to me by Prof. Larmor that an interesting 
approximate estimate of the effect to be expected might be obtained from the 
results in my paper "The Influence of Viscosity on the Oscillations of 
Superposed Fluids."f In the case of a film of oil on the surface of water, in 
which the physical circumstances are different on account of the extreme 
tenuity of the film, Prof. Lamb has made an approximate calculation of 
the effect in his Treatise on Hydrodynamics,! by supposing the surface to be 
inextensible for small wave-lengths. 

In my paper referred to, the approximations were obtained on the supposi- 
tion that the viscosity in both liquids was small — the amplitude of the waves 
being as usual very small. In the present case, although the viscosity in 
mercury is very small, that in glycerine is comparatively large. Hence a 
new approximation is now required. 

§ 2. The characteristics of the wave- motion at the free surface of a single 
liquid of great depth have been obtained for the two cases of small and of 
great viscosity respectively. In the second case there are two modes of 
vibration, both being dead-beat, for which the moduli of decay are 2hvjg and 
(O'dlvk 2 )" 1 , respectively. § For application to glycerine the. viscosity is not 
sufficiently great to make the approximation valid, except for wave-lengths 
less than 1 cm. I have thought it instructive to give this case a special 
consideration, although it is a side issue, as it affords an example of a dead- 
beat type of vibration changing to a propagated type as the wave-length is 
increased. 

* 'The Astropliys. Journ.,' vol. 29, No. 2. 

+ * London Math. Soc. Proc.,' Ser. 2, vol. 6, p. 396. 

J 3rd Edition, p. 571. 

§ Lamb, ' Hydrodynamics/ 3rd Edition, p. 567. 
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For the type of wave-motion with which we are concerned, the velocity 
system (u, v) is given by 

__ 9<£ , 9^ __ dtp ^dyfr 

ox oy oy ox 

where <f> — Ae^e tfcc+a ', 

\ 2 = Jc 2 -\-a/v ; 

v is the kinematical coefficient of viscosity and is equal to fx/p. 

Writing a\Wv = tan 2 0, the period equation can be put into the form 

tan 4 + 4tan 2 ^+^V = 4(sec^»-l)* (1) 

When v is great the term g/k d v 2 can be neglected in the determination of one 
mode, and the equation reduces to 

sec 3 + sec 2 #~f-3sec6>-~l = Q. (2) 

We obtain an admissible root sec 6 = 0*3, whence a = — 0*9 lv& 2 . 

For glycerine, v = 9*64, hence the term neglected above must be retained 
when the wave-length is greater than 1 cm. 

Writing y = sec 0, the equation (1) becomes 

y* + 2tf-4y + l+g/W>=0. (3) 

It can be solved by Descartes' method, which leads to the cubic 

K 6 + 4K 4 -4#K 2 //A 2 -16 = . (4) 

The two values of y, which are admissible, are obtained from the equation 

^-Ky + (l + iK 3 --2/K) = 0, (5) 

where K 2 is the positive root of (4), and K itself is taken positively. 

It will be seen from the table given below that the rate of decay is always 
great, but that waves of length greater than a certain limit, which lies 
between 2 and 3 cms., are propagated. 

When the wave-length is very great, one root of (4) is given by 

K 2 = — 2-f s*, where s = 4^/ Wv 2 \ 

whence y 2 — { — 2 + s%}%y + -|# = 0. 

We obtain u/k 2 v = — 2 + %i$%, 

or a = —2k 2 p±i x /gJc. 

Thus a has the same form when vie 2 is small as we obtain in the case of 
a slightly viscous liquid. 

* Cf. Basset, ' Hydrodynamics/ vol. 2, p. 312. 
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Decay of Wave-motion on 


Glycerine. 




L. 


sec 9. 


a/#*v. 


r, 


v. 
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0-3 


-0-91 


0*003 




1 


0*99 


-0-02 


-131 






i 


0-43 


-0-815 


0-013 




2 


0-89 


-0-208 


0*051 




3 




-715*0 -347* 


-0-61=1= 0-496t 


0-038 


10*0 


5 




0*92 ±0-877* 


. -0-92=1= l-61i 


0*071 


19*6 


10 




1*65 ±l-45t 


-l-40± 6'05i 


0*188 


36*7 


20 




2-95 ±3-06i 


-1 -67=1=18 -(H* 


0-629 


52*9 


30 




4-06 ±4\L8i 


-1 -96=1=38 'li 


1*18 


76*9 



L is the wave-length in cms. r the modulus of decay in seconds. 

v the velocity of propagation in cms. per second. 

§3. In the general problem we shall consider the motion of two fluids,, 
both of depth large compared with the wave-length of the motion, in the 
two cases : (1) both very viscous ; (2) the lower one slightly viscous, the 
upper very viscous. 

The period equation for waves at the interface between a liquid of 
density p and viscosity v, and a superposed liquid of density p and 
viscosity v\ is* 

4k\vp-v / P y(k^\)(k^\ , )+4k 2 ci(vp^v / p / )[p(k^x , )^p / (k^\)] 

+ p*(** + gk)(k-\') + p' 2 (* 2 -glc)(k-\) 

~ P p , [2kci 2 + * 2 {X + X r )+g'k{\--X f )] = 0, 



where 



OL 



\ 2 = k 2 +-, 

V 



V 2 



k 2 + * f . 

V 



Viscosities both Large. 

§ 4. (a) For the mode in which a is large, assuming that a is of the order 
v or v', the equation determining a is 

U^vp-v'pyik^^ik^X^^Ak^ivp-v'p^lpik^X^^p'ik^X)] 
+ [p 2 (k-\') + p ,2 (k-\)-pp' (2P + X + V)]* 2 = 0. 

The analysis is much simplified by putting v = v. With this assumption* 
after writing a — vk 2 tan 2 0, we derive the equation 

(p + p yBeo*0 + (p 2 -6pp' + p' 2 )&ec 2 0+ (3p 2 -2^ + 3/ 2 )see6^~p') 2 = 0„ 

(6) 

With this may be compared the equation (2) for a single liquid. 

* < London Math. Soc. Proc., 5 Ser. 2, vol. 6, p. 398. I have quoted it again at length 
as there are two misprints as it has been previously printed. 
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Some numerical solutions of (6) are given in the table below : — 



p7p 




0*29 
-0-91 


0-1 
0-27 
-0-93 


0'5 
0-097 
-0-98 


1 


„1 


sec 


alvTc* 





When v and v are different, we can safely say that the modulus of 
decay is not greater than (0'91v\k 2 )'~ 1 , where v\ is the smaller of the two 
coefficients. 

(b) For the mode in which a is small, we assume that a. is of the order 
v~ x . We easily find that 

- .- 9(P-P) 
2Jc(pv + p'v')' 



a = 



(7) 



This corresponds to oc = —g]2kv for a single liquid. 



Upper Liquid Viscous, Lower Slightly Viscous. 

§ 5. (a) For the mode in which a is large, writing otjk 2 v f = tan 2 6 in the 
terms of highest importance, we obtain the equation (cf 6) 

(p + p / )sec z 6+(p + p')$Qc 2 + (3p'-p)$Qoe--~(p-{-p f ) = °- ( 8 ) 

The nature of the solution of this equation is shown by the following 
table : — 



o'lo 


o-i 

0*92 
-0-15 


0*5 
0*675 

-0-544 


0*9 
0-565 
-0-68 


ajvh^ ......... 





(9) 



(fi) For the mode in which a is small we have 

2kv'p' ' 
This can be obtained from (7) by putting v = 0. 

§ 6. The Case of a Thin Layer of Glycerine superposed on Deep Mercury. — 
The period equation for this case is given in § 4 of my paper (in the first 
term of the fourth column p should be replaced by p). The effect of 
capillarity can be included by writing g-\-T2k 2 /p / instead of g in the second 
column, and g-\-Tik 2 /p instead of g in the fifth and sixth terms of the fourth 
column, where T x is the surface tension of the glycerine-mercury interface, 
.and T 2 is the surface tension of the glycerine free surface. 

For glycerine at 17°, C. v = 9*64, and for mercury, v = O'OOll, conse- 
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quently we shall approximate by putting v = at once. This greatly 
simplifies the period equation, and does not affect its validity except when 
the modulus of decay is very large. It will be further assumed that kh and 
X'h are small, so that squares of these quantities will be neglected. 
With these assumptions the period equation becomes 

(p + p' kh) u* + 4&V (p + p'kh) ** 
+ [ffp(l + kh) + F (Ti + T 2 ) + k%T 2 pjp] kot 2 
+ [^ + F(T 1 + T 2 )]4F^ 
+ k% (g + k*T 2 /p') [(p-p')ff + ^1] = 0. 

For wave-lengths less than a certain limit this equation has two negative 
roots and two complex ones. The two complex roots correspond to the 
mode in which the two liquids vibrate as a single liquid without slip at 
their common surface. The two negative roots show that the waves of the 
other mode are dead-beat. 

It appears from the numerical calculation, in the case of h = 1 mm., that 
very approximately the two negative roots are obtained from the last two 
and the first two terms of the equation respectively. Thus they are 

1 4z/[F(T 1 + T 2 )+#p] 

a 2 = — 4&V. 

These approximations are not, however, sufficiently good to be employed 
in the determination of the two remaining roots, as they necessarily make 
the real part of these roots to be positive, equal to — ^ai. 

For wave-lengths greater than the limit mentioned above all the roots are 
complex, and all the waves are propagated. This limit is given, in general, 
by the condition for equality of the roots of the last three terms of the 
equation equated to zero. 

For a glycerine-air surface, T = 65 ; for a mercury-air surface, T = 540. 
It is assumed for the purpose of an illustrative calculation that Ti = 300 at 
the common surface, while T 2 = 65. In the table below the effect is exhibited 
of a layer of glycerine of l.mm. depth on the modulus of decay and velocity 
of propagation of wave-motion at the surface of mercury. 

The damping effect of a thin layer of glycerine is very noticeable from the 
table. The rate of decay of any arbitrary disturbance will depend, in general, 
on the rates of decay of those of the modes, into which it can be analysed, 
which are propagated. The other modes disappear instantaneously. As a 
consequence of this the curvature of the upper glycerine surface is always the 
same as that of the mercury surface, a fact which may be of importance in 
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the application to an optical instrument such as the one which gave rise to 
this investigation. Prof. Wood made use of a depth of 4 mm. of glycerine ; 
a layer of this depth would give a much greater rate of decay. 



L. 


T. 


re. 


TO- 


v. 


Vo 


v . 




f -435 


0-291 










2 < 


-0026* 


-0027* 


45-5 






20-7 




[ 10 '87 


5-69 




17'5 


19-7 






-417 


0-402 










5 \ 


0*017 


-017 


303 






28'8 




L 12 *99 


10-03 




28*1 


28-9 






[ '352 


0-354 










10 < 


'081 


0-082 


1100 






39'8 




39-95 


27*03 




39-5 


40-0 




20 - 


f -518 

L t 


-518 

t 


5000 


7-19 
54-1 


7*29 
54-1 


54-1 



L is the wave-length in cms. 

t is the modulus of decay when capillarity is neglected ; r c that when capillarity is included ; 
r that of waves at the surface of clean mercury (T = 540) . 

u, v c , v Q are the velocities of propagation in the corresponding cases. 

* In the case of these quantities the approximation, on the supposition that \'h is small, breaks 
down ; hut the rate of decay is in any case large. 

f The rate of decay cannot be determined without great labour. 

If the wave-length is small compared with the depth, the rates of decay 
can be derived from (8) and (9). These approximations are applicable to 
glycerine, as the wave-length is small. The moduli of decay for the two 
modes are (l'3k 2 )" 1 and (491&) -1 ; both of these are extremely small. 

In these investigations the damping effect of the air has been neglected. 
That this is very considerable in the case of water waves has been shown in 
my paper, to which reference has been made ; but owing to the great density 
of mercury its effect will be small, especially for small wave-lengths, such as 
those tabulated. Its effect in the case of dead-beat modes is absolutely 
negligible. 

It may be mentioned that the mercury used by Prof. Wood was of com- 
paratively small depth ; this would cause considerable damping in the case of 
the longer waves. 

[Note added June 19. — It will be noticed from the foregoing table that the 
modulus of decay, when the vibrations are influenced by the glycerine layer, 
increases much more slowly with the wave-length than in the case of a clean 
mercury surface. This fact may go a little way towards an explanation of 
the damping of waves on water by oil. But a complete explanation is to be 
looked for in a totally different direction. Dr. Houstoun notes an experi- 
mental result which is apparently contradictory to usual experience.* 

* * Phil. Mag., 5 January, 1909. 
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The modulus of decay (t) for waves of length 5 cms. when the surface of 
the mercury is covered by 4 mm. of glycerine is 4*9 sees. All modes become 
propagated for wave-lengths greater than some value between 5 and 10 cms. 
For a depth of 1 mm. of glycerine this limit lies between 10 and 20 cms., and 
thus we see that as the depth is increased this critical value of the wave- 
length approaches that for waves at the surface of deep glycerine, which lies 
between 2 and 3 cms.] 



The Elastic Limits of Iron and Steel under Cyclical Variations 

of Stress. 

By Leonard Bairstow, A.R.C.Sc, Wh. Sch. 

(Communicated by Dr. Glazebrook, F.RS. Received April 20, — Read 

May 13, 1909.) 

(Abstract.) 

An explanation of fatigue is developed in the paper which is in accordance 
with all the researches on the fracture of materials by the cyclical repetition 
of stress. The theory was put forward by Bauschinger in 1886, when he 
suggested that the necessary condition of safety was that the repeated stresses 
applied should be within the limits of elasticity of the specimen, and that the 
least variation from this condition introduces fatigue and ultimately fracture 
occurs. 

In order for this to be true for Wohler's well-known experiments, the 
elastic limits must be variable within very wide limits, and the present paper 
describes observations made during the adjustment of the limits of elasticity 
to any particular condition of experiment. 

It is now found that iron or steel is capable of adjusting itself to varia- 
tions of stress, cyclically applied, after a sufficient number of repetitions. 
When the adjustment is complete, the specimen under test is found to have 
become perfectly elastic throughout the whole cycle, and fatigue does 
not occur. 

During the adjustment of the elastic limits to a given cycle of stress, 
a change of length occurs in the specimen > which corresponds to the 
extension observed in an ordinary tensile test when the yield stress is 
exceeded. For stress cyclically applied, this extension occurs even when the 
maximum stress in the cycle is less than the static yield stress. 

vol. lxxxii. — a. 2 K 



